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Abstract 

We study confinement in softly broken J\f = 2 SUSY QCD with gauge group SU(N C ) and Nf 
hypermultiplets of fundamental matter (quarks) when the Coulomb branch is lifted by small mass of 
adjoint matter. Concentrating mostly on the theory with SU (3) gauge group we discuss the J\f = 1 
vacua which arise in the weak coupling at large values of quark masses and study flux tubes and 
monopole confinement in these vacua. In particular we find the BPS strings in SU (3) gauge theory 
formed by two interacting U(l) gauge fields and two scalar fields generalizing ordinary Abrikosov- 
Nielsen-Olesen vortices. Then we focus on the SU (3) gauge theories with Nf — 4 and Nf = 5 flavors 
with equal masses. In these theories there are Af = 1 vacua with restored SU (2) gauge subgroup 
in quantum theory since SU(2) subsectors are not asymptotically free. We show that although the 
confinement in these theories is due to Abelian flux tubes the multiplicity of meson spectrum is the 
same as expected in a theory with non- Abelian confinement. 
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1 Introduction 



According to Mandelstam, Polyakov and 't Hooft jjj confinement of charges arises as Meissner effect upon 
condensation of monopoles. Once monopoles condense the electric flux is confined within the electric flux tube 
P connecting the heavy trial charge and anti-charge. The flux tube has constant energy per unit length 
- the string tension T = (27ra') _1 , and this ensures that confining potential between heavy charge and anti- 
charge increases linearly with their separation. However, since dynamics of monopoles is hard to control in 
non-supersymmetric gauge theories, this picture of confinement for many years remained to be an unjustified 
qualitative scheme. 

The breakthrough in this direction was made by Seiberg and Witten in 0, |^. Constructing exact solution 
to Af = 2 supersymmetric gauge theory they have shown that the condensation of monopoles really occurs near 
the monopole point in the moduli space of the theory once Af = 2 supersymmetry is broken down to TV = 1 
by the mass term of adjoint matter ||]. After the work of Seiberg and Witten it has become very important 
to understand to what extent this Abelian confinement of electric charges is similar to confinement of color 
we expect (but cannot control) in real QCD. Moreover, one expects the QCD-like confinement also in Af = 1 
supersymmetric QCD which can be obtained as a limit of large mass of the adjoint matter \i of softly broken 
Af = 2 QCD, but again, we have no control on the dynamics of this theory in the large limit. 

One important distinction noticed by Douglas and Shenker || appears in SU (N c ) gauge theories with 
N c > 3. Since SU(N C ) gauge group is broken down to [/(l)^^ 1 by the VEV's of adjoint scalars there are 
iV c — 1 generically different flux tubes, one per each U(l) factor. Numerous flux tubes lead to existence of 
too many hadronic states in the spectrum |(| (see also (T) for the D-brane reinterpretation of this result). In 
particular, there are N c different sets of quark-antiquark meson Regge trajectories Q For example, in pure 
SU (N c ) gauge theory the number of families of sets of these trajectories with different slope is the integer part 
of (N c + l)/2. The presence of many quark-antiquark meson trajectories reflects the essentially Abelian nature 
of confinement in Seiberg- Witten theory. 

On the other hand in Af = 1 supersymmetric theories described in the framework of Seiberg's duality Q 
we should get non- Abelian confinement since without adjoint fields breaking of the gauge symmetry down to 
Abelian subgroup does not occur. It is believed that condensation of "magnetic quarks" of non- Abelian dual 
theory leads to confinement of ordinary quarks. Still the mechanism of this confinement is not understood. The 
problem is that usually non-Abelian dual gauge groups, say SU(r) groups do not admit flux tubes. 

The bridge between two approaches was suggested in (9[ |l0| . It was shown that certain Af = 1 vacua of softly 
broken Af — 2 QCD at small \x preserve non-Abelian SU(r) subgroups of SU(N) gauge symmetry (r < N c ). 
The theories at some of these vacua correspond at large [i to Af = 1 QCD described by Seiberg's duality. This 
suggests that confinement at these vacua should be non-Abelian even at small /i. In particular we should have 
only one set of quark-antiquark meson Regge trajectories. In this paper we are going to study the outlined 
above proposal for the non-Abelian confinement from the side of softly broken Af = 2 QCD at small [i. In 
particular we analyse what happens to Abelian flux tubes when a non-Abelian subgroup of gauge symmetry 
is restored at certain Af = 1 vacua. We will mostly focus on theory with the SU(3) gauge group and Nf < 5 
flavors of fundamental hypermultiplets to be called as quarks below. In sect. || we review the vacuum structure 
of the theory studied partially in |j| px[ . We consider the case of large quark masses tua, A = 1, . . -Nf, to 
keep theory at weak coupling and discuss vacua where quarks develop VEV's. These vacua are classified by the 
number r of different colors for which quarks have non-zero VEV's. For SU(3) gauge theory the isolated vacua 
can have r = 0, 1, 2 and we are interested mostly in weak coupling vacua with r = 1 and/or r = 2 for large 
values of tua- 

In sect. |3| we present weak coupling Abelian description of the low energy effective theory around these 
vacua and discuss the low energy spectra. Generically the SU(3) gauge group is broken down to its maximal 
U (1) x U(l) Abelian subgroup by the VEV's of adjoint matter and we have two light photon multiplets. In sect. [| 
and sect, g we study in detail magnetic flux tubes in r — 1 and r = 2 vacua responsible for the confinement 
of monopoles. It turns out that r = 1 vacua possess the standard Abrikosov-Nielsen-Olesen (ANO) flux tubes 
§ | described by one [7(1) gauge field and one scalar field. In the limit of small masses of the adjoint matter 
[i these flux tubes are BPS-saturated [jlT) . However, at r = 2 vacua the flux tubes turn out to have more 

1 Each set corresponds to main trajectory together with "daughter" trajectories. 
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complicated structure. They are formed now by two gauge fields interacting with two scalar fields. We consider 
the lattice of these flux tubes of more general type and study which ones among them are BPS states. To do 
this we calculate the interaction potential of strings at large distances. 

In sect. U we consider finally the r = 2 vacua in the SU(3) gauge theories with Nf — 4 and Nf = 5 in 
special limit of equal quark masses, when the SU{2) subgroup of gauge symmetry is restored. The reason 
is that classically restored SU(2) subgroup stays unbroken also at quantum level since corresponding SU(2) 
subsectors are not asymptotically free (though, of course, the SU (3) gauge theory with Nf = 4 and Nf = 5 
is asymptotically free itself). We study what happens to the Abelian flux tubes and monopole confinement 
under restoration of the SU (2) gauge subgroup and find that one of flux tubes becomes unstable and eventually 
disappears from the spectrum. Two other "elementary" strings become essentially identical in this limit. This 
eliminates unwanted multiplicity of the hadron spectrum mentioned above since in this vacuum we get only one 
set of meson Regge trajectories as expected in a theory with non-Abelian confinement. 



2 Vacuum structure 

2.1 Superpotential and vacuum equations 

Consider softly broken J\f = 2 SU(N C ) gauge theory with Nf flavors. The field content in terms of M = 1 
supermultiplets can be described by the gauge multiplet W a and chiral multiplct $ in the adjoint representation 
(forming together J\f = 2 gauge supermultiplet) so that first contains the gauge field {A^- and complex Weyl 

fermion (A a )*., while the second - complex scalar and complex Weyl fermion {ipa)], = 1, ■ ■ ■ ,N C ), - 
all being N c x N c matrices with zero trace and 2Nf chiral multiplets of matter Q A ' % and Qa,% in the N c and 
N c representations respectively, i.e. i = 1, . . . , N c and A = 1, . . . , Nf. The vacuum structure we are going to 
discuss can be associated with extrema SW of the superpotential 

Nf N a -l 

A=l ' k=2 

Here the last term generally breaks J\f = 2 supersymmery down to J\f = 1 and we have chosen it to be linear in 
the independent invariant combinations Tr$ fc . In fact, later on we simplify it further assuming that only the 
coefficient fi2 = M m front of the mass term of adjoint matter is non-zero. 

To find vacua of the theory we also have to impose the D-term condition which reduces to [<&\ $] = and 

Nf / -At \ 

»■ X Q A,i QAj - Q ' Qa,3 = v# s (2) 

A=X ^ ' 

The parameter v = v{Q,Q,Q,Q) = jj-TrD is in fact some function of quark's VEV's and it can be nonzero 
only if all Qj ^ for j = 1, . . . , N c . The variation of the superpotential ([!]) gives rise to the following set of 
equations 

V2$>)Q A - J +m A Q A ' 1 -0 VA = l,...,N f 

(3) 

V2QA,i^j + m A QAj = VA = l,...,N f 
including also F-term condition J^- = 0, which should be taken into account together with (g). One should also 
remember that Tr$ = 0, say, introducing the Lagrange multiplier /ii into (^), the same is to impose vanishing 
condition only onto the traceless part 



Sit _„,- \\ ^ . / ,.l nt 5 





f - 1 E Q Ak Q^ +E fc ^ - ) = w 



of V2F* + E fc > 2 fc^fc (t> k ~% with 



Nf 

M = J] Q A -Q X j . (5) 



/ j ^1A,J- 
A=l 
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Let us first investigate the structure of solutions to vacuum equations (0), (||) and We shall basically use 
quasiclassical regime when all masses of the quarks tua S> ^qcd where the coupling is small and quasiclassics 
is a good approximation. Whenever it is not possible we will add extra arguments based on the Seiberg-Witten 
analysis of strong coupling regime of SUSY gauge theories [Q, || . 

Equation [$^,$] = gives rise immediately to the conclusion that, in general position, up to gauge trans- 
formations 

/ fa . . . \ 











V o ... <t> Ne J 

Tr$ = J2 4H = 



(G) 



and gauge group SU(N C ) is broken down to U(l) Nc 1 since the off-diagonal vector fields {A^. acquire the 



masses proportional to 



[$,A4 = (^-^)( J 4 /i );. 



(7) 



due to the Higgs effect. When fa = fa for some i ^ j non-Abelian gauge symmetry may be partially restored. 
For diagonal $ @ the second term in (Q) vanishes for i ^ j and one ends up with a simple problem of counting 
of the eigenvectors of matrix We shall consider in detail the case of gauge group SU(3). 



2.2 SU(3) gauge group 



The SU(3) Af = 1 SUSY gauge theory without matter has exactly N c = 3 vacua, all in the strong coupling 
regime. This is a particular case of general situation for an SU(N C ) pure gauge theory, which has N c points in 
the moduli space when N c — 1 monopoles become massless (see, for example, 0, |l3[ ||]). 

Let us add one flavor with mass m,A — rn so that the equations (||) and linear combinations of the equations 
(|J) turn into 

(V20i + m)Q 1 = Qi{V2fa + m) = 



(V20 2 +m)Q 2 = O Q 2 (V2fa+m) = 
-V2fa ~ V2<p 2 + m)Q 3 = Q 3 (-V2fa - V2fa 2 + m) = 



\/2(oiQ 1 -g 3 Q : 

V2 (q 2 Q 2 - Q 3 Q :: 



(8) 



-(2 M 2-3M302)(201+0 2 ) 

-(2 A12 -3m3<?!>i)(0i +2</> 2 ) 



where we have restricted ourselves to (the only essential in SU(3) case) nonzero /x 2 and ^3. Matrix $ may have 
only one eigenvector, which can be by gauge transformation turned to any given direction in colour space, for 
example n 



and 



Q 



A,i 



Q % = uS\ = u 





\0 J 



(9) 



Qi = udj = u ( 1 ) (10) 

Taking into account (JsJ) , ( ^0| ) equations (||) give rise to \f2fa +m = or fa = — -^=. Then, from the last equation 
of (||) and vanishing of Q 2 and Q3 one gets that fa = — ^- = and, as a consequence fa = —fa — fa — j^=. 

2 In the case of the SU (3) gauge group we will use common notations u, d and s for the Q l with i = 1, 2, 3. 



4 



As a result matrix <!> acquires the form 



$= — 

V2 



—m 











rn 





T 






rn 














~2 



(11) 



It means that the gauge group is in fact broken only up to St/ (2) x t/(l). The vector fields A 12 and A* 3 (and 
their superpartners) become very heavy, of the mass ~ m 3> A, while the mass of the photon A 1 ^ and its 
superpartners is totally defined by "higgsing" via Qi and, from the fourth equation of d8t), is given by 



light 



flu m fi 2 



Ay/2 



fi 3 m 



(12) 



i.e. is of the scale of SUSY breaking for generic \ii and ^3 Q. In ( |12|) we have also introduced the coupling- 
constant g which comes from kinetic term and corresponds to correct normalization of mass as a pole in the 
propagator. From the equations (@) and formulas (||) and (|l0|) we immediately find the value of the quark 
condensate (QQ) in this vacuum 



m I u 2 =u-tm 

2 V 4v^ 

(dd) = (is) = 0, 



c 

2' 



(13) 



where we introduced parameter £ = 2(uu) The SU(2) subgroup is unbroken at the scale of order of m, it 
interacts only with very heavy matter Q2 and Q3 (of m A) and by the Seiberg-Witten mechanism it runs to 
the strong-coupled phase and has two Seiberg-Witten vacua j|, || . It means that vacua of the full theory would 
correspond to 

/ -m \ 



V2 



V 











^±A 



SU{2) 











m 



(14) 



where A^ (2) m = 



sc/(2) y 

= A 5 . The gauge fields A 2 ^ and their superpartners would get masses of the order 

of Asu(2) "C rn and the unbroken t/(l) subgroup of this SU(2) will remain light, with the mass ~ y/ H2A.su {2)- 
This mass is even less than (|l2|). In order to avoid running of the SU(2) subgroup into the strong coupling one 
has to add more flavors to our gauge theory. 



2.3 More flavors with SU(3) 

Let us now turn to the situation with more flavors. Of course we can still have vacua of the type considered 
in previous section for each flavor, however, in this case more complicated vacua when different flavors get 
simultaneously non-zero VEV's also arise. 

From (|3|) it follows that either Q A and Qa are zero for a given A or they are eigenvectors of the diagonal 
matrix (^|) with eigenvalue —niA/V^- An important restriction on extra vacua also comes from the formula (|J) 
with i=£ j (see @) 

N f 

/■;; X Q A,( l\; Mi (15) 

A=l 

3 In this paper we do not consider possible complications related to the "fine-tuning" of (^), say when for [13 = 4 ^^ 2 extra 
fields could become massless, at least classically. 



Combining these conditions together we see that for example, in the case Nf = 2 one may consider only 

qA,i = Qi § Ai^ j e 

/ 1 \ 



Q 2,« = d 2 p 2 = d 2 





V o J 

(°\ 

1 

V o / 



(16) 



and Qj 



Qi$Ai, or 



and equations (pi) give rise to V2~4>i 



Qi.i = UiSn — u\ ( 1 ) 
Q 2 ,r = d 2 S i2 = d 2 (0 1 ) 
mi = and \[2fa + m 2 = implying 
/ mi \ 



(17) 



1 

71 



m 2 
\ 





-mi - m 2 J 



(18) 



This is a general situation up to gauge rotations in the colour space, so one gets an extra vacuum with two U (1) 
gauge groups softly broken at "light" level of (|l2|). Indeed, from (16), (|l7j ) and, following from (||) relations one 
finds the following values for the vacuum condensates: 



1 3 
(uiu 1 ) = ^=(2^2 - 3^ 3 2 )(2<?!>i + fa) = (fj, 2 + — -j=li3m 2 )(2mi + m 2 ) 



V2' 



(19) 



{d 2 d 2 ) = t=(2^2 - 3^ 3 (^i)(0i + 2fa) = (^ 2 + — 7=^3Wi)(mi + 2m 2 ) 



V2 



2V2' 



The spectrum of light fields will be discussed in detail in sect. |3.3| . The total number of vacua is in agreement 
with the formula from [fl0[ 

mia(No-l,Nf) 

#(vacua)= (N c -r)C? f (20) 

r=0 

where r counts the number of nontrivial eigenvectors or solutions to (^|). Indeed, r = term gives N c (= 3 
for SU(3)) Seiberg-Witten vacua "without matter", r = 1 term adds {N c ~ 1) • N f (= 2N f for SU(3)) vacua, 
corresponding to (|IJ) for each flavor) - with the gauge group broken to SU (N c — 1) at the scale m ^ A in weak 
coupling. The term with r = 2 



{N c - 2)C. 



N f (N f - 1) 



(21) 



AT C =3 2 

corresponds exactly to the situation, considered in this section. 

Formula (po|) has a simple physical meaning. The factor (N c — r) is exactly the Witten index of unbroken 
by the adjoint matter gauge group (in the case we consider adjoint matter always breaks SU (N c ) down to some 
"lower" SU(N C — r)). The combinatorial factor C r f counts the number of possibilities to arrange quark VEV's 
within Nf flavors. 



2.4 Colliding vacua and Higgs branches 

If some of the masses of the matter multiplets coincide (ttia = hb for A ^ B; A, B = 1, . . . , Nf) one gets a 
Higgs branch with VEV's (Q) ^ 0, "growing" from the corresponding point on the moduli space of the Coulomb 
branch. The (real) dimension of the Higgs branch is 47i, i.e. the dimension of the hyper Kahler manifold is 
four times the "number of hypermultiplets" TL. 
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In the r = 1 case one can always choose the nontrivial eigenvectors of $ along some fixed direction in the 
colour space, say only Q 1A ^ and Qai ^ for any A. Then we get 4Nj real parameters (for coinciding all 
jtia = m and matrix <J> of the form of (pi])), which should obey two real F-term and one real D-term relations. 
One extra degree is "eaten up" by the 17(1) gauge group so that finally the dimension of the Higgs branch 
appears to be AN f - 3 - 1 = 4(JV> - 1) or H r =i =N f -l. 

In the r = 2 case the situation is a bit more complicated. The solutions for nontrivial Q and Q can be now 
chosen in the form of the following rectangle matrices 



and 



\Q 



Ak\ 



( Q 11 


Q 12 ... 


Q 1Nl \ 


Q 21 


Q 22 ... 


Q 2N f 


\ 





... 









( c 


?11 Ql2 


\ 






c 


?21 Q22 


















(22) 



(23) 



V QNfl Qn s 2 / 



consisting of the eigenvectors of $ — diag(— to, — m, 2m), with 8Nf parameters. These parameters obey now 
four real D-term conditions (|J) and eight real F-term conditions and gauge group SU(2) x U(l) "eats" four 
"phases" , so that the total number of independent parameters is 8Nf — 12 — 4 = 8(Nf — 2) or H I= 2 = 2(iV/ — 2), 
which coincides with the formula of H for r = 2. 



2.5 Baryonic branches 

Formula ( pp| ) is based on counting of the number r of nonzero eigenvectors - solutions to (^|) . We have considered 
the cases r = 1 and r = 2 for the SU(3) gauge theory, but it is obvious that the maximal number of eigenvectors 
is r = N c = 3. This possibility appears first time for Nf = N c (=3 if we still consider the SU(3) gauge group) 
and is only possible, however, in the situation when Xm=i mA = Ea=i m A = 0. 

Consider, for example, r = 2 vacuum of the previous section. The relation ([l5]) still requires Q A,% = Q l 5 Al 
and Qa,i = Qi^Ai, but now for i, A = 1, ... ,3. In this way we get the set of vacua, parameterized by 



(qa,iQ A1 -q a , 3 q a ' 

1=1,2 

£ (q a ,2Q a < 2 - q a , 3 q a - 3 



'' ] -- uiu 1 - s 3 s 3 



(M2 



A=l,2 



= d 2 d l - s 3 s = (fi 



2V2 
3 

2V2 



A*3?7i2)(2mi -I- to 2 ) 
/x 3 toi)(toi -I- 2to 2 ) 



(24) 



The system (Q) describes the baryonic branch ^|, indeed (|T^) is simply obtained from (24) putting Q3 = Q3 = 0, 
and in this point TO3 can be "unfrozen" from —(mi +TO2). This is a baryonic branch since the VEV's of baryons 



B = u'd 2 * 3 = ^e ljk e ABC Q A ' 1 Q Bj Q C ' k 
B = u x d 2 h = ^£ ljk £ ABC QA,iQB,jQc,k 



(26) 



are nonzero. The dimension of this baryonic branch is #(Qi, Qj) — #(i ? ) — #(£ ) ) — #(phases) = 12 — 4 — 2 — 2 = 4. 



To "symmetrize" (E4J) one may also add a relation 



{Qa,iQ A ' 1 - Qa,2Q A ' 2 ) = uiu 1 - d^d 2 = (/X2 H ^pM3"i3)(mi - m 2 ) 

A=l,2 

which is not independent, but is just the difference of two in (p^j). 



(25) 
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The simplest analog of this branch exists already for the SU(2) gauge theory with Nf = N c = 2 matter 
hypermultiplets with mi = —mi ~ m. The corresponding moduli space is described by 



(27) 




and 



Q feA = 



QAk = 




(28) 



or eight real parameters and the relations 
(the last one 



2|) and (Q) give one and two real relations on them correspondingly 
Q 1 Q 1 - Q 2 Q 2 = 2^m (29) 



is a particular case of counted carefully general situation when the i — j relations of (Q) impose 2 • rank = 
2 • (N c — 1) conditions) plus a U (1) phase, so that the dimension is 8 — 2— 1 — 1 = 4. The corresponding baryon 
operators are B = Q 1 Q 2 and B = QiQ 2 - 

One may also add more flavors, up to Nf — 4 (the conformal point for the SU (2) theory) when the baryonic 
branch can be described in terms of (|27]) and 



Q 



kA 




!Ak 



(30) 



Q 3 

V o g 4 / 

modulo (||) and (||) which gives 16 — 2 — 1 — 1 = 12 parameters or three hypermultiplets. 



3 Low energy mass spectrum 

In this section we work out the action of the effective low energy Abelian theory and use it to study the low 
energy spectrum in isolated Af — 1 charge vacua of the theory at vha ^> A (i.e. at different values of quark 
masses when there are no Higgs branches). We put jUfc = Li5k2 or keep only fi 2 = A* 0; assuming also that 
[i < A. 



3.1 Abelian description 



Let us consider scales of order y/JImX which are well below W-boson masses at small fi. At this scale SU (3) 
gauge group is broken down to U(l) 2 by VEV of the adjoint scalar (0) 



<l> = 



f 01 \ 

02 

V o o 03 / 



o 





-a 3 



Q8 

V3 



-2 



«8 



= A 3 a 3 + A 8 a 8 



(31) 



V3 ) 
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Figure 1: Root and weight vectors in the Cartan plane for SU(3) group. We have depicted explicitly the root 
vectors oi and 02 (the simple roots), the highest root ol\2 = 0:1+0:2 and the weights of the fundamental 
representation \i x = \x u , fi 2 = an d A*3 = Ms; corresponding to u, d and s quarks respectively. We have also 
depicted the "normalized" roots (of unit length) ei = 012/V2 and e2 = OL2/V2 as well as normalized weights 
u, d and s of the length 1 / a/3- 



at generic values at quark masses. We also have two U(l) gauge fields (in the orthogonal basis denoted A 
and A^'; note that the orthogonal basis is normalized so that Tr<f> 2 = h (a 2 + <z|)) introduced via 



(3) 



/ A {8) 



A, 



\ 



A 



(8) 



V3 



.4 



(8) 



V 



a 3 4 3) 



As4 8) 



(32) 



V3 / 



where our notations correspond to expanding gauge and adjoint fields in the orthogonal basis of the Gell-Mann 
matrices - the basis of diagonal Cartan generators of the SU(3) Lie algebra (see (149) and (15C) in Appendix). 
It will become clear later that it is more adequate to use the basis given by (divided by two) root vectors in the 
Cartan subalgebra of SU(3) 012/2 and 02/2 [| see Appendix and fig. |l], but we find it technically simpler to 
work in the orthogonal basis. The price to pay for this is non-integer electric charges of u and d quarks with 



respect to orthogonal gauge fields, say A^ and A)?' , see below. 

In these notations the bosonic part of the low energy effective Abelian theory acquires the form 



,(8) 



Sqed 



d A x 



1 



¥ 2 



1 



4.9 2 



1 2 1 

-? |<9 M a 3 | + -x |<V 
9 9 



N f 

E 



E 

1 i—u.d., 



(33) 



+ V(u, d, s, a 3 ,a 8 ) 

5 We will also used "normalized" roots eo = <y.\j\f ! i e2 = cc2/v / 2 and ei = CX12/V2 and "normalized" weights, e.g. u = pj/x/2. 
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Here u, d, s label different colors of scalar components of quark supermultiplets Q l while 



\/3 v/3 p 



(34) 



where we introduced also (non-orthogonal) components of the gauge fields interacting directly with u- and 
d-quarks 

4«) = ^Uw + 



2 ^ 



(35) 



The potential in the Lagrangian 



f E 



|£«l 2 



<9W 



0*r 



V(u, d, s, a) 



5 2 D 



<9W 



(9a 



2 ^/ 



<9W 



dQ A ^ 



aw 



i<fc 



E E 

A=l i=u,d,s 

2 — dW f dW 

2 t^UW 



<9a 3 



aw 



<9a 8 



E E 



i—u,d,s 



dW 



dW 



dW 



dQ 



A,i 



(36) 



is given by D and F terms (g), (|)) (see also formulas (||), (f|)), expressed in terms of the scalar components 
Q 1 = u, Q 2 = d, Q 3 = s and a 3y8 (3|) 



v(u, d, s, a 3 ,a 8 ) = ?-Yl "j"! • ff 2 



- D d d f + |! (d: + - 2/^) 2 



' 2 



F" - + y/2, 



Nf / 
A=l \ 



a 3 H — + V2mA 
v3 



2 , 3 2 



-L (f« + F d d - 2f;) + V2 ( 



-a 3 H — + V2m,A 
v3 



(l< 



A I 2 



-^a 8 + V^m^ 
v3 



(l^ + lSxl 2 ) 



Id 



(37) 



where color indices i,j as in (|33j ) run over the set of u, d and s. 

The QED coupling constant g in ( [33] ) is small at small /i, determined by the mass of light quarks and photons 
(cf. eq. ©) 

(38) 



1 fim 



With logarithmic accuracy we do not distinguish between two different coupling constants associated with two 
U(l) factors in (|||). Below we consider different types of N — 1 vacua and analyse what kind of flux tubes do 
we have in each type. 
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3.2 r=l vacua 



Consider first r = 1 vacuum when only one Q A ' 1 = Q 1 ' 1 quark flavor develops VEV. We drop flavor index below 
in this subsection because all other r = 1 vacua have the similar structure. As we have learnt in sect. 2.2 the 
VEV's of scalar fields in this vacuum are given by ( |l3| ) where we assume for simplicity that \i and m are real 
and positive and choose the phase of the it-quark condensate to vanish. The adjoint scalar develop VEV given 
by (11), in terms of fields as and as this reads 



(«3) 



2a/2 
2V2 7 



(39) 



From the kinetic term for u-quark in ( |33| ) one can get the mass matrix for the U(l) gauge fields (A^ , ) 



M~ = -g urn 



1 



1 1 

\ 4 3 



whose eigenvalues are 



together with 



= 2g 2 fj, m = -g 2 C, 



mi=0 



(40) 

(41) 
(42) 



where we have used parameter £ = 2(uu) — 3/im. These equations show that one massive photon appears since 
it-quark develops nonzero VEV jl3| ) and breaks one of two U(l) gauge groups while the ma ssles s photon is 
associated with the second unbroken U(l) group. This fact we have already pointed out in sect. |2.2|. Classically 
the SU (2) group which includes the latter U(l) factor remains unbroken, see for example formula (|ll|). However, 
in quantum theory this SU(2) subgroup is broken down due to the Seiberg-Witten mechanism |5) and the 
second photon acquires small mass of the order of yJnKsu(2) due to monopole/dyon condensation in this SU (2) 
sector. 

Now let us consider the mass matrix for the scalars and as, to the leading order in /i/m it can be read 
off the two last lines in potential ([37]) ■ In this approximation it coincides with the photon mass matrix ( f40| ) 
and therefore one complex a-field remains massless while the other one acquires the same mass ( |4l] ) as one 
of the photons. This is directly related to the fact that N = 2 supersymmetry is not broken in the effective 
low energy QED ([33j) in the leading order in fi/m ^ [l4j. To see this note, that in this approximation the 
perturbation of superpotential (|I|) proportional to fj, is linear in fluctuations of a-fields. Thus it boils down to 
the Fayet-Iliopoulos (FI) F-term which does not break M = 2 supersymmetry in the effective QED Q . In 
the next to leading order in fi/m (which corresponds to taking into account fluctuations of the a-ficlds in the 
third line of potential j37|)) TV = 2 supersymmetry is broken and Af — 2 supermultiplets split 14 . Below we 



restrict ourselves to the leading order in /i/m so that Af = 2 supersymmetry is preserved in the effective low 
energy Abelian theory (|33|). 

To complete the study of the mass spectrum in vicinity of r = 1 vacuum consider the mass matrix for 
quarks. The d- and s-quarks are very heavy (with masses of the order of m > A) and, hence, decouple from the 
low-energy theory (see last two terms in (p7|)). The mass matrix for remaining four real components of u-quark 
(Reu, Rcu, Imu, Imu) can be read off the terms 



a 2 a 2 
6 1 u> 2 



F: + y/2, 



^a 3 



9_ 
2 



+ %/2/jas 

v3 




i (^) 2 
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of the potential (j37j) and has the form 



1 



M 2 u = 2g 2 fim 



\ 



10 

o o i i 



0-- 



(43) 



with one zero eigenvalue corresponding to the state "eaten" by Higgs mechanism, and three other eigenvalues 
coinciding with the photon mass (^) . These three real states of it-quark combine with two real states of massive 
a-field and three states of massive photon to form the bosonic part of one long Af = 2 supermultiplet which 
contains eight boson and eight fcrmion states. The long multiplet appears because electric charges are screened 
in the broken U(l) sector by Higgs condensate, therefore the corresponding central charges of Af — 2 algebra 
vanish, and short BPS multiplets cannot appear p| . To sum up, one gets one long Af = 2 multiplet with mass 
( pl| ) formed by the u-quark condensation and another one with much smaller mass of the order of ~ \/ ^su(2) 
formed by the monopole/dyon condensation in the classically unbroken SU(2) sector. 



3.3 r=2 vacua in the low energy effective theory 



Let us consider now the r = 2 vacua of sect. 2.3. Assume for simplicity that we have only two quark flavors 
with masses mi and m 2 . In this vacuum quark fields develop VEV's ([19]) while adjoint VEV's are given by 
([l8|). In 03, as basis they are 

(03) = 7=( m l ~ m 2) 



V2 



(44) 



(mi + m 2 ) 



The mass matrix for the gauge fields (A^ , ) can be read off the kinetic terms for u- and d-quarks in ( |33"1) 
and has the form 



Mi 



2 



1 + UJ 



1 



1 



V V3 



^(1 + w) 



where we have introduced the parameters of r = 2 vacua 

£ = 2(uiu 1 ) = 2//(2mi + m 2 ) 

and 



(d 2 d 2 



2too 



mi 



{UlU 1 - 



2mi + rrii 



analogous to the parameter £ of r = 1 vacuum ( J13| ) . Two eigenvalues of this mass matrix are given by 

K) ± - g -Y^ ± 



with 



n± = - (l + uj ± ^fl-uj+u) 2 



(45) 

(46) 
(47) 

(48) 
(49) 



We see that for generic values of mi and m 2 both U(l) groups are broken and both photons acquire masses. 
The mass matrix for two complex fields 0,3 and ag, is identical to ( p!5"| ) in the leading order in \ijm as can be 
seen again from (|3~^). 

The mass matrix for quarks is now of the size 8x8 including four (real) components of w 1 -quark and 
four components of d 2 -quark. It has two zero eigenvalues associated with the two states "eaten" by the Higgs 
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mechanism for two U(l) gauge factors and two non-zero eigenvalues coinciding with photon masses (Hq). Each 
of these non-zero eigenvalues corresponds to three quark eigenvectors. Altogether we have two long Af = 2 
multiplets with masses (ff8|), each one containing eight bosonic and eight fermionic states. 

Now let us briefly comment on the special case of coinciding masses mi = (cj = 1) to be considered 
in detail in sect. |6[ In this case (03) = and SU{2) subgroup of the SU(3) gauge group is restored on the 
Coulomb branch at zero \i at least classically, see (|18|). However when one switches on the terms proportional 
to \i this SU(2) subgroup becomes broken completely by u- and <i-quark condensates. It is easy to see that 
all three masses of SU(2) gauge fields are the same and given by ([f8| ) with fl + = 2 which corresponds to the 
special value u> = 1 in ( |49|) when mass matrix becomes diagonal. 



4 Flux tubes 

In this section we consider flux tubes in isolated Af = 1 charge vacua of the theory at ^> A and \i <C A. In 
this regime the low energy effective theory reduces to Abelian model (^) at weak coupling so one can use the 
semi-classical methods to study string solutions. 



4.1 ANO strings in r=l vacua 

At r = 1 vacuum d- and s-quarks are heavy and we can simply ignore them. The VEV's of light fields are given 
by eqs. ([l3|) and (^), in particular u- and it-quarks have the same VEV's. This suggest that one can look for 
the ANO string solution using the ansatz (|l3|), (|39|), i.e. to fix 

d=s=d=s=0 



3 y/3 ( 5 °) 

— =777, a 

>V2 

and express u and u in terms of a single complex field tp 



a-i = =772, ax = =777 



^/2 



(51) 



Since u-quark interacts only with particular combination of the gauge fields and A^ given by ( ^5| ) 
it is natural to assume that only this combination is non-zero on corresponding string solution. To implement 
this let us rotate the fields Ajj, and A^ to another orthogonal basis 



,(«) _ y^3 4 ( 3 ) 1,(8) 

a* 2 ^ 2 tJl 

A 2 „ = --A^ + 

2 2 



(52) 



The meaning of the subscript "2" of means that it is directed along 0.2 ~ e.i (see fig. |l|) and at the moment 
we only need that 

A^ = (53) 



on our string solution. It is easy to check that the ansatz (|5_y); (|53| ) satisfies the equations of motion. With this 
ansatz the bosonic part of softly broken Af — 2 QED ([53]) reduces to the form of standard Abelian Higgs model 
(the relativistic version of the Landau-Ginzburg model) 



*ah = / d 4 , (f^) 2 + |vH^| 2 + A (M 2 - c) 1 



(54) 



with particular value of quartic coupling and electric charge n e — 1/V3 of the field tp, see (|34|). The field <p 
develops VEV dl|) |<^)| 2 = C = 3|/jm|, therefore U(l) gauge group is broken, or the photon acquires mass 
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m 2 = 2g%/3 while the Higgs mass is equal to m H = 4A£. Note, that the photon mass coincides with that of 

([ll]) , since the fields A^/ and A2 ll diagonalize the photon mass matrix (^0|) , if it were not so the fields A^f* and 
^2^i would mix leading to contradiction between condition ([53]) and equations of motion. 

Strictly speaking the substitution (03) = — ^/^ m an d ( a s) — ~^^ m does not satisfy equations of motion 
following from (|33|). In fact VEV's 03 and a§ ( |50| ) get x-dependent corrections of the order of y/ fx/m H. 



However, as we already explained in sect. 3.2 one can neglect this effect in the leading order in /i/m and consider 
£ to be just a constant £ = 3/im. In this approximation the perturbation term in the superpotential ( p]) i s linear 
in a and reduces to the FI i^-term which does not break M = 2 supersymmetry in effective QED W, |14 ], As is 



explained in detail in |L4j the ANO strings in U(l) theory are BPS-saturated in this limit (see also ||19[). Below 
we use the term "BPS string" for these "almost BPS" solutions, they belong to the short JV = 2 multiplets 
which become long J\f = 1 multiplets when next to leading order corrections (breaking Af = 2 supersymmetry) 
are taken into account |14]. 



Let us now briefly remind the basic features of the BPS ANO strings [jri|. For arbitrary A in (|34|) the Higgs 
mass run (the inverse correlation length) and photon mass m 7 (the inverse penetration depth) are different and 
their ratio is important parameter in theory of superconductivity, characterizing the type of superconductor. 
Namely, for mn < rn y one has type I superconductor, while for ran > to 7 it is of type II, this is related to 
the fact that scalar field produces an attraction for two vortices, while the electromagnetic field produces a 
repulsion. The boundary separating superconductors of the I and II type corresponds to ran — m~ l , i.e. to 
special value of quartic coupling A 

x ^ntf = £ (55) 
2 6 

This is exactly the value of quartic coupling in ( pi) ) we get from the potential ( |37| ) using the ansatz (|50|), (|5l|). 
In this case vortices do not interact. It is well known that vanishing of interaction when vtih = m-y can be 
explained by the BPS nature of the ANO strings. The ANO string satisfies the first order equations and saturate 



the Bogomolny bound 11 , which can be found from the following representation of string tension T (see (|54|)), 



(M 2 -C) zu) 2 + \\" 



T = 2ni\n\+ / d 2 x J (J-Fu ± ^= (M 2 - C) e u ) +^\Y lr -±i, u V., r -\ 2 } CXil 



Here indices I,J= 1,2 denote coordinates transverse to the axis of the vortex. For positive n we take the upper 
sign in (|56|), whereas for negative n we take the lower sign. The minimal value of the tension is reached when 



both positive terms in the integrand of (56) vanish, then the string tension becomes 

T BPS = 27rC|n| (57) 

where the winding number n counts the magnetic flux 2im. The linear dependence of string tensions on n is 
consistent with the absence of string interactions. 

For the n = 1 case vanishing of the integrand in (p&) leads to well-known two first order differential equations 



r^(r)-/(r)0(r) = 



(58) 



r ar 3 

(for the positive signs in (pq)), where the profile functions 4>(r) and f(r) are introduced in a standard way, i.e. 

4 u \ x ) = V3e u ^ [f(r)-l] ^ 

Here r = yjj2j=i 2 X J ^ s tne distance from position of the vortex and 1? (pi'd — —eijxj/r 2 ) is polar angle in 
transverse to the axis of vortex (1, 2)-plane. The profile functions are real and satisfy the boundary conditions 

0(0) = , /(0) = 1 

(60) 

</>(oo) = VC> /(oo)=0 
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which ensures that scalar field reaches its VEV \/( at infinity and vortex carries one unit of magnetic flux. 
Equations J5E| ) with boundary conditions ( |60| ) lead to unique solution for the profile functions (although an 
analytic form of this solution is not found, though for £ = the system (^) is equivalent to the "radial" 
Liouville equation). The tension of string with winding number n = 1 is given by particular case of ([57]) 



T n= i = 2ir£ = 67t|/ito| 



(61) 



In Af = 2 QED emergence of the first order equations ( pq) means that some (half) of the SUSY charges of 
N = 2 algebra act trivially onto the ANO solution (cf. Jl6|, |17|, |Tg| , Q). In this case the Bogomolny (topological) 
bound for the string tension coincides with the central charge of SUSY algebra. 



Now let us discuss the embedding of flux (winding) numbers of the string (p0|), (51), ( p9| ) into the Cartan 
subalgebra of SU (3) group. Throughout this paper we will use the convention of labeling the flux of a given 
ANO string by magnetic charge of the monopole which produces this flux and can be attached to its end 0, this 
is possible since both string fluxes and monopole charges are elements of the group 7r 1 (t/(l)® 2 ) = Z® 2 . This 
convention is convenient because specifying the flux of given string we automatically fix the charge of monopole 
confined by this string. 

9|) has non-zero while ^2 M vanishes, i.e. the matrix for this 



The string solution (50), 
configuration looks like 



A, 



A 



V3 











o -- 



0-- 



(62) 



It preserves the SU(2) subgroup acting in the 2x2 right low corner of this matrix. Thus the flux of this string 
(charge of the monopole attached to its end) is orthogonal to the root vector a 2 ( |152[ ) (or e2, see fig. [j]). Hence, 
the string charge vector is proportional to u-quark weight vector on the Cartan plane, fi u ~ fi 1 . We call this 
string u-string and for this reason denote its charge vector q u . 

To work out the absolute value of the it-string charge we recall that the Dirac quantization condition in the 
r = 1 vacuum with non-zero VEV of the u-quark looks like 



n 
2 



Thus for winding number n — 1 one has 



(63) 



(64) 



where we use that |u| = n e = l/VS, see (|34|). 

From (|6l|) we see that absolute value of the string charge is |q u | = -s/3/2 and this vector points into the 
middle of the side of adjoint hexagon, see fig. ^. Note, that monopole charges (whose positions in the Cartan 
plane coincide with the positions of W-boson charges) correspond to the corners of adjoint hexagon. We see 
that the component of monopole charge parallel to vector u is confined by the ANO u-string. Moreover, it is 
easy to see from fig. ^ that the values of projections of charges of ei and eo monopoles onto the u-direction 
exactly matches the charge of the u-string. Thus u-string confines u-component of ei and eo monopole charges. 
This remark can be used as an explanation of rather strange from group theory point of view value of charge 
of the u-string lying outside the root lattice of SU(3) group. 

However, all monopoles present in the theory (ei, e2 and eo- monopoles) have also non-zero component in 
the orthogonal direction to vector u, i.e. along e2. What happens to this component of the monopole charge? 
The answer to this question can be found in quantum theory, where the story is a bit more complicated. 
As we already explained in sect. 2fl the SU(2) gauge subgroup is preserved only classically at r = 1 vacua, 
and from the Seiberg-Witten exact solution [Q, |5| we know that in quantum theory SU (2)-subsector runs into 
strong coupling regime where the SU(2) subgroup is never restored; instead, one gets either monopole or dyon 
vacuum. Consider, say, the monopole one. At nonzero /i monopole develops VEV of the order of nAsu(2) and 

B This monopole can be a superposition of "really existing" monopoles in given vacuum as we will see below. 
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Figure 2: Strings and different phases in r = 1 vacuum. Black squares correspond to the monopole and quark 
states in confinement phase while white squares - to the states in Higgs phase. Vectors q M and q2 = e 2 /2 label 
the (charges of) magnetic and electric strings correspondingly. 



this monopole has charge e 2 , see ([L4[). It means that e 2 -components of monopole charges are screened by the 
condensation of the e2-monopole. 

Passing to dual magnetic theory at strong coupling one can also study the formation of the ANO electric 
string due to condensation of e2-monopole along same lines as above. It is clear from the Dirac quantization 
condition that e2-string arises with the charge q 2 = ^2/2, see fig. 0, It confines e2-components of d- and 
s-quark electric charges while their u-component is screened by the w-quark condensation. Note again, that 
e2-components of d- and s-quark electric charges exactly coincides with the charge of the e2-string. The electric 
e 2 -string which confines the fundamental charge is very similar to strings at monopole vacua of SU(3) gauge 
theory without fundamental matter (labeled by r = in our notations) studied in fl5| . 

To sum up, in each r = 1 vacuum we obtain one BPS magnetic ANO u-string and one BPS electric G2- 
string. In other words we have mixed Higgs phase for u-quark and e2-monopole and confining phase for ei- and 
eo-monopole as well as for d- and s-quark, see fig. 0. 



4.2 Strings in r=2 vacua 



Now let us turn to strings in r = 2 vacua. For simplicity we assume, first, the presence of only two quark flavors 
with masses m\ and 1112 under the same conditions as in sect. B.3, i.e. we consider s-quark to be heavy and 

and (44), in particular VEV's of u- 



ignore it in low-energy theory. The VEV's of all fields are given by (fL9| 
u-quarks are equal to each other as well as VEV's of d- and d-quarks 
using the following ansatz 

u =u A = d A ,i—^ 



and 

Therefore, we look for string solutions 



d A =d A = 5 A :. 



V2 



(65) 



while the fields 03 and a§ in the leading order in /x/m are given by their VEV's. 
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With this ansatz the effective action (133) becomes 



£Vl 2 - 



0d \ 2 - e (1 - c)) 2 + £ (i^i 2 + i^i 2 - e(i + ^)) 2 ) 



(66) 



4 d^i 2 



where we use the same notations as in ( p4[ ) and (35) while £ = |<y5 u | 2 and cj = | c/?^ | 2 / 1 | 2 are introduced in ( ^6| ) 
~~7|) (note that we consider both of them real and positive). Gauge fields A^jf 1 and Ay here as functions 



and 



of orthogonal fields Ay and Ajf J are given by (|35|) . To study more general case of complex £ and u> one should 
modify ansatz (|6^) taking into account relative phases of the fields Q and Q. 

Let us now derive the Bogomolny bound for string solutions in the theory (|66|). Assuming that all fields 
depend only on two spatial coordinates orthogonal to the string axis one can rewrite ( |66|) as follows 



1(8) 



T = I d 2 x 



+ 



2<? 



p(8) , 5 



1 

+ 2 



n( u ) _L_ ■ V7( M ) 



(|^ 1 | 2 + |^| 2 -£(l+cj))e /J + 

(u,d) 



(67) 



where we introduced (two-dimensional) dual field strength for gauge fields A)? 1 ' as F^ u,d ^ = ^eijFjj' ' . The 
upper sign here corresponds to positive total flux, the value of J d 2 x + F^u^j, while the lower sign 

corresponds to the negative total flux. 



The last term in fl67| ) represents exactly the string fluxes while all other positive terms in (|67]) should be 
zero on the BPS solutions, leading to the following first order equations 



1 



F, 



(3) 



I J 



2ff 



-F, 



(8) 



4V3 



e(i^i 2 -i^i 2 -e(i-^))e/./ = o 

e(\vu\ 2 + \w\ 2 -Z(l+u J ))e I j = 



(68) 



7 («) 



ieeijVy ip u = 



V^ d) + ieeuVy ifd 







where e = ± is the sign of total flux. 

One can classify possible solutions to these equations by behavior of the fields at spatial infinity in (1, 2) 
plane. One type of solutions has nontrivial winding for the u-quark another type has winding for the ci-quark 
and there are also mixed solutions when both u- and d-quarks wind at infinity. 

Let us start with string solutions when only u-quark winds. Assuming for simplicity the unit winding number 
one has at x — > 00 

(69) 



This behavior requires the following behavior at infinity for the gauge fields 



1 



V3 1 



and 



4' 



(70) 
(71) 
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which ensures finite contribution to string tension coming from kinetic terms of u- and d-quarks. Such behavior 



at infinity means that the flux of the field A^ vanishes while the flux of the A^' field is 

V3. 



2tt 

for unit winding number, then representation (|67]) for the BPS string tension gives 

Ti = 2tt£ = 4?r \n{2m x +m 2 )| 

To implement conditions (Fll) let us rotate the gauge fields A^ , A^ to another orthogonal combinations A 



(72) 



(73) 



(given by (|35|)) and Ai^ (the component along the vector ei) as 

= ^f, 3) + ^ 



.4 



4 d) 



2" -" 
V3 



A (3) 

2 



2 " 
f -AW 



(74) 



and rewrite our first order eqs 



The behavior at infinity (|6 
functions 



in terms of these fields. 
7l| ) suggests that one can look for a solution in terms of the following profile 

<Pu(x) = 0„(r)e 1 '', ifd(x) = 4> d {r) 
A u (x) = 2ejj ^ [Mr) - 1], Af\x) = VSejj ^ f d (r) 



(75) 



Note, that the charge of M-quark with respect to Ai M field is n e = 1/2, while the charge of d-quark with respect 
1^ is n e = 1/V3, which exactly corresponds to the coefficients in (|7^). With this substitution the first 
order equations (|6^) turn into the system of four first-order nonlinear differential equations 

r^Mr)-fi(r)Mr) +\fd{r) <f> u {r) = 



dr 



(76) 



r dr 



3 



Mr) 2 -\Mr) 2 + \t-<"t) = <> 



The reason why the fields — an d A}? with trivial behavior at infinity cannot be simply put to zero is 
mixing between ip u , A\^ and tfd, A^ since the fields (p u and Ai M which wind at infinity are not eigenvectors of 
the mass matrix ( p5[ ) . In other words ( |76| ) does not have nontrivial solutions with fd = and <pd = a/^? (which 
is a naive ansatz suggested by the boundary conditions (69), (fn])). 

First order equations (|7^) should be supplemented by the boundary conditions which ensure that quark 
fields go to their VEV's at infinity, u-quark field has no singularity at r — (4> u (0) — 0) and A\^ field carries 
flux equal to 2vr (/ x (0) = 1). 

Now let us discuss the charges of our string. Since A^ goes to zero at infinity the gauge field matrix looks 

like 

/ 1 \ 



A„ 







(77) 



Al£ 

r— *oo Z 

\0 -1J 

at infinity. Eq. J7^ ) means that the charge of string is directed along the root a 12 ~ ei in the Cartan algebra. 
Its charge vector should satisfy the Dirac quantization condition qiu = 1/2 (u-quark winds) which gives 



Qi = ei 



(78) 
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Figure 3: Lattice of string solutions in r = 2 vacuum. We have specified explicitly (1,0) = ei, (0, 1) = e 2 , 
(1, —1) = eo, (1, 1) = ei + e 2 and (2, —1) = 2u strings. 



This is the reason why we call this string e!-string. Its charge is directed to the upper right corner of the adjoint 
hexagon, see fig. H. 

In a similar way one can consider the string solution when d-quark winds at infinity, while u-quark field runs 
to its VEV. This gives us e2-string with charge vector directed to the upper left corner of the adjoint hexagon, 
see fig. [| The easiest way to describe this string is to use another set of orthogonal gauge potentials A 2fi and 
related to A lfi and A*jf* by 

A (u) = ^_ Al _ l A (d) 

* (79) 
A 2 » = \a x , + ^4«> 

The behavior at infinity of e2-string is given by 

fd(x) = <j) d (r)e ld , y u {x) = <Pu(r) 
A 2I (x) = 2e u % (f 2 {r) - 1) , A™ >(x) = VSe u ^ / u (r) 



The profile functions here satisfy equations and boundary conditions similar to ( [761) . The tension of the e2-string 
is given by 

T 2 = 2ttujZ = 4tt \fj,(2m 2 + m x )| (81) 

and this string, if exists, is also BPS saturated. 

In general, one can also consider string solutions when both u- and d-quarks wind. Let us fix some basis, 
say, the potentials Ai^ and Affl d72) and introduce the profile functions for the string when u-quark winds 
n times and d-quark winds k times (to be denoted as (n, fc)-string since its charge in our normalization is 
qn,fc = nei + ke 2 ) defined as 

<Pu{x) = K{r)e m \ = Mr)e M , 

A u {x) = 2ee u ^ (A(r) - A(0)) , Af{x) = yfi se u ^ (f d (r) - f d (0)) 
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These profile functions satisfy the first order equations (J76|) , the only modification is due to the extra sign factor 
e in (|8^) which accounts for possible positive or negative total flux 

n + ujk 

£ = £ n ,k = i ; n = sign(n + Luk) = ±1 (83) 

\n + u)k\ 

The boundary conditions for profile functions of (n, fc)-string read 

MO) = 0, <M0) = 0, u (oo) = ^/e, <M°°) = yfct (84) 

for the scalar fields and 

for the gauge fields. The boundary condition for /i(0) is easy to derive from / u (0) = e n .k n using relations (|79|). 

From the boundary conditions ( |85| ) and two first equations of ( |76| ) one immediately finds the leading behaivior 
of the profile functions at r — > 

<j) u (r) oc r fl ~ fd/2 = r e "' fc ™ 

~° , (86) 
<j>d{r) oc r fd = r e "' fcfc 

which means that nonsingular at the origin solutions to the BPS first-order equation exist only for both n and 
k either positive or negative simultaneously (dependently on the sign of t n ,k)- 

The tension of (n, fc)-string is given again by last term in the representation ( |67| ) and reads 

T n ,k = 2<|n + ujk\ = 4?r \h {n{2m 1 + m 2 ) + k(2m 2 + mi))| (87) 

The lattice of possible BPS (n, fc)-strings is shown in fig. ||; ei-string is identical to (1, 0)-string while e2-string 
is nothing but (0, l)-string. There is one more string whose charge belongs to the adjoint hexagon, namely the 
(1, — l)-string (to be called also e -string in what follows), see fig. [|. As is clear from fig. ||, ei-string can be 
considered as a bound state of eo and e2-string, which exactly corresponds to decomposition of the root c*i2 
into linear combination of the simple roots a.\ and a^, see fig. [y. It is clear, thus, that string charges are labeled 
by (normalized to unity) root vectors of the SU(3) root lattice. We will return to this question in sect. [|, where 
we study which of the BPS (n, fc)-strings exist as solutions to the first order equations ( f76|) and which among 
them are stable. Here we conclude by considering certain special cases in which some of (n, fc)-strings becomes 
ANO strings and equations ( |76|) reduce to ANO equations (|5Sl ) with one gauge and one scalar potential. 

Consider, first, the simplest example of this kind with both quarks having unit winding number, namely the 
(1, 1)- or (ei + e2)-string. This string can be considered as bound state of ei- and e2-strings. If this bound 
state exists and saturates the Bogomolny bound then (K7]) suggests for its tension 



T Xt i = 127r \^(m 1 +m 2 )\ 

Now let us give some evidence that this BPS string indeed exists. Consider the special value of parameter ( f47| ) 
u> = 1, then VEV's of u- and d-quarks are equal. The string solution can be easily constructed with the only 
component of gauge field being non-zero (putting — 0) and one complex scalar field introduced via 

<Pu =<Pd=^j= (89) 

With this ansatz the theory (|6^) reduces to standard Abelian Higgs model of type (|54|) with equal values for 
photon and scalar masses m 7 = run = <? 2 £/3 and this model obviously possesses the BPS ANO string solution. 
This mass of the Abelian Higgs model coincides with the eigenvalue of the photon mass matrix, corresponding 
to f2 = 2/3, see ([48|), (|49|). This is one of possible explanations why a single gauge potential and single scalar 

(8) \ k 

field appear in the string solution at to — 1 : gauge field ^ and scalar ( |89[ ) correspond to eigenvectors of gauge 
and scalar mass matrices. 
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One can reach the same conclusion directly from equations (f76|). Substituting the ansatz c/) u = (f>d and 
fx = 3/2 fd (this is consistent with boundary conditions ( p5| ) because /i(0) = 3/2 and /d(0) = 1 for (1,1)- 
string) we see that four equations ( ff^ ) at w = 1 reduce exactly to the system of two ANO equations (|58|), under 
the following modifications: the r = 1 parameter £ is replaced by the r = 2 parameter £ and the coefficient in 
front of the second term in the last equation of ( |5§| ) replaced by m^/2£ = g 2 /6 for a> = 1. 

Now let us move parameter u> away from u> = 1. It is clear that under continuous deformation string solution 
cannot disappear or become non-BPS state, since BPS string belongs to a short multiplet which cannot turn 
into long one without breaking of some amount of supersymmetry (the number of states cannot jump). Thus, 
at least at some re gion of values of parameter u> we expect existence of the BPS (ei + e2)-string, we specify 



this region in sect. 5.2. where more detailed analysis of the nature of (n, fc)-strings and their stability, based 



on study of string interactions, will be presented. One may also use more direct evidence based on numerical 
simulations (we have done this using the MAPLE program) of the BPS first order equations @ Q 

To conclude this section note, that electric flux tubes quite similar to our magnetic (n, /c)-strings were found 
in Af = 2 SU (3) theory without fundamental matter at strong coupling vacua |Q . These strings appear due 
to the condensation of monopoles/dyons and their charges should be associated with the (normalized) weight 
lattice on the Cartan plane instead of (normalized) root lattice which arises for our magnetic strings. Strings 
of this type arc often called Z3 strings when studied directly in the underlying non-Abelian ST/ (3) theory, see 

M Eil El PI M I2I. 



4.3 BPS formula for the string tensions 



In this section we rewrite the mass formulas for the BPS strings studied in previous sections in a form which is 
more familiar for the BPS objects. Tensions of BPS strings are given by central charges of J\f = 2 SUSY algebra 
@ 0, H, [h§. They have general structure 



T BPS = 2n\qJ\ 



(90) 



where, in our particular case, q s is the charge vector of a given string in the Cartan plane (see fig. ^ and fig. ||) 
and f is the generalized vector parameter of the FI _F-term, which can be defined as 



where (see @), (pjft and (152)) 



f = -Afi(f> 



4> = 4>\OLx2 + 4>20<-2 = V2(0iei + (/> 2 e 2 ) 



(91) 



(92) 



where <pi and 4>2 are the corresponding components of the VEV's of the adjoint scalar matrix (^) in given 
vacuum. In the orthogonal basis of fields 03 and a§ they are equal to (see (|3"T|)) 



1 



«3 



-0,3 



08_ 

V3 



(93) 



Let us now show that (|90|) indeed gives correct tensions for all BPS strings considered before. Start with r = 1 

(94) 



vacua, from ([TT|) one finds that 



= _ m ( ei -|) 



In r = 1 vacuum we have only u-string with the charge |u, see (J64]). For this particular string ( |90|) gives 



T u = 8vr 



e 2 

— jiui u I e! 



= 67r|/im| 



(95) 



cf> = — miei — m 2 e 2 (96) 

7 We tried to use the string solution in monopolc vacuum presented in Jlit , but unfortunately it does not seem to satisfy all four 



which coincides with (|6l|). Now consider r = 2 vacuum, where one has (see 



cqs.(M) 
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thus, for the (n, fc)-string formula ( |90| ) gives 

T n ,k = 87r|^i(nei + fce 2 )(miei + m 2 e 2 )| = 47r |/x (n(2mi + m 2 ) + k(2m 2 + mi))| (97) 

and this result coincides with (|87j). The BPS formula ( |90| ) is valid in the limit /i — > 0, £ ~ /im^ = const, when 
one can neglect breaking of N = 2 SUSY in effective QED ([33|) and it assumes also the weak coupling regime 
\m A \ > A. 

We conclude this section noting that masses of W-bosons are also determined by the BPS formula, following 
from (0) 

m w = V2\cx W (p\ = 2\qw<j>\ (98) 

where qw is the (normalized) charge of corresponding W-boson in the Cartan plane. In particular, it means 
that tensions of eo-, ei- and e2-strings which belong to the adjoint hexagon are proportional to the masses of 
corresponding W-bosons at least at weak coupling when mw ^> A. Indeed, charge vectors q s of the eo-, ei- 
and e2-strings coincide with the charges qw of related W-bosons, and comparing (^) with (^8|) we see that 
tensions of all strings are proportional to the masses of corresponding W-bosons. 



5 (n,k)-strings 

In this section we address the question of existence and stability of (n, A;)-strings in r = 2 vacua studying their 
interactions at large distances. This allows us to see which strings attract each other and form stable bound 
states and which do not interact so that their bound state is only marginally stable. First, we consider the 
interaction of the ANO strings in order to develop necessary technique and then turn to the interactions of 
(n, fc)-strings. 



5.1 Interactions of ANO strings 

Let us first develop the method of effective vertex to calculate the interactions of ANO strings. This method is 



well-known in instanton physics and is used there to calculate instanton interactions |26| , |27j , 28 , |29|. Here we 
generalize it to the case of solitonic ANO strings (it also can be easily generalized to any solitonic branes). 

Transform the scalar and gauge fields of the ANO string from the "regular" into "singular" gauge making 
the U(l) gauge transformation with exp (— i{f), where # is polar angle in the (1,2) plane orthogonal to the 
string. This gauge transformation is singular at r — so now the topological charge of the string (flux) comes 
from small circle around the origin instead of the large circle at r — * oo, and in this gauge the substitution ( p9| ) 
turns into 

tp(x) = 4>(r) 
A i{x) = —eu — f(r) 

where the dependence on arbitrary electric charge n e is restored. Here profile functions <j)(r) and f(r) satisfy 
the ANO first order equations (cf. with (|5^)) 

m - f{r)cj>{r) = 
ar (100) 

r ar 



and boundary conditions 



were n is integer winding number, while 



0(0) = 0, /(0) = n 

(101) 

i(oo) = /(oo) = 



e n = 



sign(n) (102) 



n 
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Using equations (100) which determine the exponential fall-off of the functions / and <fi ~ V? a t infinity, 
corresponding to the Yukawa behavior in two transverse dimensions, we get for the large r asymptotic 



Ax(x) 



1 



Xj _ 

e/j —y/m^r e 



(103) 



where c„ is the coefficient to be fixed below, while e n is given by ( |102| ), both c n and e„ depend on the winding 
number n. We will also need the behavior of the (tw o-dimensional) dual field strength F = \tijFu at the 
infinity in two-dimensional plane (1,2), and from (103) one finds 



F 



(104) 



Now let us work out the effective vertex for the ANO string. This vertex once added to the tree level QED 
action (e.g. (p4[ )) should reproduce all effects of presence of the ANO string in the framework of perturbation 
theory (see ]26|, B71 Eq, p~~ 



For the string with winding number n we propose the following form ^| (up to the 



overall normalization factor) 

V^ NO = J L>A(<x)expj- J d 2 a2y/ 7 hic n 



<p<p(X(tr)) 



2n e g 



;F^(I(o-))n M1/ 



(105) 



where fields ip and F^ are considered as functions of string co-ordinates X{a). We also introduced here 
an antisymmetric tensor n liV (X{tr)) orthogonal to the string world-sheet at the point X(cr). For the case 
of the straight string at rest directed along the third axis n M „ = for fj, or v = 0, 3, while nu — eu for 
fi,i> = I, J = 1, 2. Note, that here we use static parameterization for the string world-sheet, a% = t and &i = x^, 
in other words 

dX a dX p 



1 

2 £lJ da i da j 



(106) 



and this interaction is nothing but well-known interaction of string with antisymmetric tensor i?-field B = *F. 
To check the expression for the effective vertex ( jl05| ) let us calculate the correlation function 

■/string 

(107) 

in the string background assuming that all points x's and j/'s are far from the axis of the string at X{<r). On 
one hand this correlation function can be calculated in scmiclassical approximation just substituting classical 
expressions (103) for the scalar and gauge fields at large distances from the string axis into the correlation 
function (107), on the other hand the same correlation function can be calculated in perturbation theory as 

r ANO 



^{x)...F{y)...V.r U ) 



(108) 



once the effective vertex V^ NO is added to the tree level QED action, and this should give rise to the same 
result as a substitution of classical fields. 

To see this expand scal ar fi eld around its VEV ip = i/f + Sip and write down the bilinear in scalar fields 
term in the exponential in ( |l05| ) as 



ipp 



ip- 



(109) 



Consider, first, the linear in quantum fluctuations terms. Expanding the exponential in ( |108| ) in Sip and in F we 
calculate the correlation function (108) to the leading order in coupling constant using the tree level propagator 



(5ip(3 



d 2 a5p{X{a))) - — K {m 7 r) 

Z7T 



2V27T y/f 



(110) 



8 In what follows we are going to keep only bosonic background fields in our expressions for the string interaction, the fermionic 
terms can be restored by supersymmetry, similar to the instantonic case, see for example [ESl . 
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for the massive scalar and 



(F(x) / <faF(X(a))) 



9 



d 2 



2ir dr 2 



K (m^r) 



g 2 m 2 e ~ m -,r 



2V2tt 



(111) 



for the massive vector field, where r is the distance between point x and position of string in (1, 2) plane. These 
correlation functions are nothing but propagators of two-dimensional theory rewritten in four dimensional 
notations. 



The qua dratic in quantum fluctuations term in (109) cannot be verified using just tree level approximation 



(110), (111), since it is next to leading effect in coupling constant. Still one can restore the quadratic dependence 
on field ip in the exponential in ( |l05|) observing that full effective vertex can depend only on fields (p and cp 
rather than upon their VEV's. The result looks a bit surprising from the point of view of string theory, and 
the origin of this quadratic dependence is that in our picture string tension T = {2ira')^ 1 is "dynamical" and 
determined, in contrast to the fundamental string theory, by condensate of a scalar field. 

Let us now determine the constant c„. Ignoring quantum fluctuations in (105) and substituting the scalar 
field by its VEV v^? one gets 2-\/27rc n ^ in the exponential in (105), which should be equal to string tension 
27r£|n| (cf. @); that gives 



fmi 



Substituting this into (|105| ) we finally obtain 



V r f NO = J DX(<t) cxp j-2^M J d 2 a 



2n e g 2 



and the effective partition function of low-energy QED is now given by 



J DA^DipDipexp I ~S QED - ^V r 



ANO 



(112) 



(113) 



(114) 



The effective vertex in the exponent allows to have many strings with different winding numbers, hence theory 
in (114) should be treated perturbatively since non-perturbative effects (strings) are already taken into account 
explicitly. 



Let us now use the effective vertex (113) to compute the interaction potential of two straight static strings 
directed along the third axis at large separation R in (1, 2) plane. Expanding exp (— V^ N °^ in powers of 
Vn NO we keep only the term V r ^ 1 NO V^ s NO , related directly to the interaction potential of strings Ui t % via 



(yANOyANO\ 
\ n\ ri2 I 



j DX 1 {a)DX 2 {(T)e- 



U!, 2 V2 



(115) 



where Xi(cr) and X2(tr) correspond to two strings while V% is the volume of the two dimensional space in 
(0, 3)-plane. Using the propagators (110) and (111) to calculate the correlation function in the l.h.s. of (115) 
we finally obtain 

C/i.2 = -(27r)gflni||nal . [1 - £i£ 2 ] (116) 

where E\, £2 refers to the signs of winding numbers ( |102| ) of two strings. 

We see that potential of string interactions has exponential fall-off at large separations R. The first term in 
square brackets comes from the exchange by scalar field, and it always gives the attractive contribution to the 
potential. The second term comes from the photon exchange and its sign depends on the relative signs of ri\ 
and U2- If ^1^2 > this term gives repulsion which cancels attraction produced by scalar exchange, which is, of 
course, a well known result - the BPS ANO strings do not interact. If, however, n\ri2 < (which corresponds 
to string-antistring interactions) the photon exchange also gives an attraction. Thus the total string-antistring 
interaction potential is always attractive. 
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Let us conclude this section presenting the effective vertex of the ANO string which takes into account its 
motion, ft is clear that such generalization of (|113|) has the form 



V, 



ANO 



J £X(<T)exp(-S,f in s) 



string = 27r | 



V det g ind ipip(X(<r)) 



2n e g 



r F pl/ (X(cr))n A , I/ 



Here gi n d is the determinant of the induced metric given by 



ind 

9u 



dX^ dX^ 



(117) 



(118) 



dai da j 

and interaction with i?-field is defined in (|106| ). Besides the usual N amb u-Goto term the string action contains 
also higher derivative terms omitted in (117). The effective action (117) takes into account interactions of the 
bulk fields and <p with the two-dimensional "field" X fJi (cr) living on string world sheet. 



5.2 Interactions of (n,k)-strings 



In this section we apply the effective vertex method to calculate the interaction potential of (n, fc)-strings in 
r = 2 vacua. First, from the first order equations ( |7r^ ) we get the behavior of the fields <j> u (r), 4>d{r) and fi(r), 
fd{r) at infinity 



M r ) = \/I( 1-c n, 



Mr) 



0j£ 1 + Cn,k 



Fx{r) 



3(0 — |w) y/myr 



(119) 



Fa(r) 



2m^ 



where c n .k will be determined below, while sign factors e n> k are given by 
gauge and expressions for dual field strengths F\, Fd related to fi, fd via 

£n,k df(r) 



n e r dr 



In ( |119| ) we use the singular 
(120) 



F 

where n e = \ and n e = for /i and fd respectively. 

Clearly the leading asymptotic behavior in (119) is determined by the lightest gauge and scalar fields so the 
photon mass in (HE) coincides with smaller eigenvalue in formulas (PH), (49), i.e. m 1 = (m 7 )_ and CI = £!_. 
In fact ( |119| ) already determines the relation between the fields (p u , fd, Fi and Fd and eigenvectors of the mass 
matrix (g5|) ip^, tp(+\ F^ and 

6ip u = cos (3 v 6 (p(~~> + sin P 9 8ip( + ' 
Sipd = — smPipSip^' + cos (3 v 8ip( + ' 



(121) 



and 



(+) 



Fi = cos PfF^ + sin (3 F F 
F d = - sin fipF^ + cos f3 F F (+) 



(122) 



since (3 V and (3f are determined by ratios of the coefficients in front of leading asymptotics in (HE) 

+ a 2 

Ian U v - - — - 

tan/Jp = 



Cl 



(123) 



Vsci 
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Of course, the same relations can be derived directly by diagonalization of the corresponding mass matrices, 
see sect. B.3. Using these formulas and asymptotic behavior in (119) for the lightest scalar and vector fields at 
large r, one gets 



i{3 v ^fm 



F ( \r) = e n , k c n ,k 



2m 2 



1' 

ra^v 



(124) 



^COS/3f y 7 ^ 



+ 



It is clear that the effective vertex for (n, fc)-string would effectively depend in the leading order only upon the 
lightest scalar and gauge fields and ' since only these fields determine the large r be havior of string 
interaction. Follo wing t he sa me steps leading to effective vertex for the ANO string in sect. 5.1 and using again 
the propagators (110), (111) we arrive to the following expression 



V n , k ~ J DX(<r)exp^- J d 2 o2^c n , k 



dip 



(-) 



e jhh. F (-) n 



(125) 



up to normalization factor. Let us finally fix the constants c n ,k here. To do this rewrite expression (^7|) for the 
string tension in the form 



T n ,k(<f) = 27r \n\ip u \ +k\ip d Y 



(126) 



Expanding in (126) the scalar fields around their VEV's one can extract the linear term in 5<p( ' using relations 
(121). Comparing it with the expression linear in ^i^ - ' in the exponent in (|125|) we find for the coefficients c n ^ 



'' j3 v (n - fctan/^v^) 



(127) 



Substituting this back to (125) and taking into account ( |126| ) we finally arrive to the following effective (n, fc)- 
string vertex 

\/det3i„ d |n|(^ u | 2 + k\ip d \ 2 \ 



DX{cr) exp 

cos 2 (3, 



-2ir 



d 2 a 



g 2 flcosf3F 



(n - k tan f3 v \fui) F^Jn^ 



(128) 



where, as in ( 117 ) we have already taken into account motion of string. Here only F^ u enter the exponential 
because only this component survives at large distances. However, on general grounds it is natural to expect 
that the effective vertex depend on the following combination of two gauge fields: q ni fcF where we defined 
vector field strength as 

Fp, = eiF llu , + dy/3F$. (129) 
Taking this into account we can conjecture that the effective vertex has the form 



V n , k = / DX(tr) exp \ -2tt / d 2 a y/dctg md \n\ip u (X (<r))| 2 + k\^ d (X((r)f 



1 cos 2 (3, 



g 2 Q cos 2 (3 f 



P q n ,fc • F fiU (X(a))n flI/ 



(130) 



Exponential here depends only on the component of the field strength, directed along the charge of the 
(n, fc)-string. The extra piece that appears in (13C) as compared to (128) associated with component of 
the gauge field a nd d o not contribute at large distances. We are not going to use the effective vertex in the 
(conjected) form (130) below and write it down only in the sake of completeness. 

Let us finally compute the interaction potential of strings at large distances, technically for this purpose it is 
easier to use effective vertex in the form (125). Expanding exp ^— J2 n k Vn,kJ in powers of V n .k we keep again 
only the term (V ni _k 1 V n2i k 2 ) ■ Calculation of this correlation function to the leading order in coupling constant 
using tree level propagators (11C), (111) leads to the following interaction potential of (n, fc)-strings 

lu) 2 + \uj 2 e -™-rR 



16 — (Q 



£l£2 



(131) 
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where c\ = ct ni ,ki) f° r (ni, fa)-string and c 2 = C(n 2 .k 2 ) f° r ( n 2, fa)-string are given by (127). Quite similar to the 
case of ANO strings the first term in square brackets here comes from the exchange by two scalars tp u and (p<i and 
this interaction always gives rise to attractive potential. The second term arises due to exchange by two photons 
Ai^ and the sign of this contribution is determined by product of the sign factors ( |S3| ) for two strings. 

In particular, for [n\ + ujki)(n 2 + ujk 2 ) > the photon exchange gives repulsion which cancels the attraction 
produced by the scalar exchange and in this case strings do not interact. If instead (ni + ujk±)(n2 + ujk 2 ) < 
the interaction potential is attractive so that two strings form a stable bound state. 



5.3 Lattice of (n,k)-strings 



Now let us use the interaction potential (131) to verify the existence and stability of the (n, fc)-string solutions. 
We are going to determine which strings among the solutions to the BPS first order equations (f76|) are stable 
and which do not exist as BPS states. The similar approach was used in |50|, [u| to study BPS soliton states in 
two dimensions and various dyon states in Seiberg-Witten theory. 

Suppose first, one has two BPS strings (ni,fa) and (ri2,fa) with (m +uiki){n2 + u>k 2 ) > 0. Then according 
to ( |131| ) these strings do not interact, it means that the bound state of these strings (m + 7i2,fa + fa) exists 
as a BPS solution to the first order equations (f76|) but it is only marginally stable. The tension of this string is 



given according to (B7J) by sum of the tensions of its " components" 

T le2 = T 1 +T 2 (132) 

where T s = T(n,,fe s ) []■ H instead (m + u>k\)[n 2 + wfa) < 0, two components attract each other and form a 
stable bound state. However we do not know a priori if this bound state is BPS saturated or not, all we know 
is that its tension is within the bounds 

|Ti-T 2 | <T im <T 1 +T 2 (133) 



where the lower bound is given by the BPS formula (p7|) while the upper bound follows from the potential (131) 



From (|33|) we see that interaction of strings depends on the value of parameter uj. Therefore let us first 
distinguish physically different regions for this parameter. The masses of W-bosons are classically given by ( |98| ) 
and, as follows from (0), ( fTcf ) and ([l7]), they are proportional to 

(134) 





M 






Mz\, 




1 


M { Z -\) 




z 







where instead of mi and m.2 we have introduced new variables M and z defined by 

M 2 = (2mi + m 2 )(2m 2 + mi) 



(135) 

z = UJ 



One can see from (134) that at special values u> — 0, uj — 1 and ui — 00 one of the W-boson masses vanishes and 
it corresponds to the restoration of corresponding SU (2) gauge subgroup. However, as we mentioned before this 
is correct only classically and in quantum theory the SU (2) subgroups are never restored for the theories with 
Nf < 4 (in the next section we consider the theories with four and five flavors where the SU(2) subgroup is 
restored at u> = 1). Instead the SU (2)-subsector runs into strong coupling and W-bosons never become massless 
[|[ ||. Hence, it is clear that three regions around values w = 0, 1 and 00 are in fact strongly coupled and we 
cannot use our semiclassical analysis there, therefore one has two separated weak coupling regions 

< uj < 1 and 1 < uj < 00 (136) 

and we cannot pass from one region to another within weak coupling regime keeping real uj. However, in 
the effective low-energy Abelian theory these three points look differently: at uj = and uj = 00 one of the 

9 Note that we restrict ourselves only to real and positive values of parameter u> (^t]). For complex ui one would have far more 
complicated lattice of [n, fc)-strings. In particular the mass formula ( pTy suggests that we probably can have stable BPS bound 
states at complex u) which then decay at real u). Thus the surface Imo; = should be the curve of marginal stability (CMS) for 
these bound states. 
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Figure 4: Lattice of string solutions for < u> < 1. The non-BPS strings are drawn by white dots while the BPS 
strings are depicted by black dots. All BPS strings are marginally stable bound states of "fundamental" ei- 
and e2-strings, corresponding to two "closest" roots of SU(3) algebra. The non-BPS strings can be "crossed" 
by straight line n + cuk — or e(w) when one moves parameter u> in the region < u < oo. 



photons becomes massless (see (JIg|)), which is not true for oj = 1. In fact, in theory with classically restored 
SU{2) subgroup value u = (u> — oo) corresponds to the Argyres-Douglas point |52| [33| |34|, |3^] where the 
monopole/dyon vacuum (r = 1 vacuum) collides with charge vacuum (r = 2). Thus the values u> = and 
to = oo correspond to strong coupling even in the Abelian low energy theory. On the contrary, nothing special 
happens in the Abelian theory ( p3|) at u = 1, the low energy theory does not "feel" the restoration of ST/ (2) 
subgroup. This means that in the low energy description we have only one region < to < oo instead of two 



regions ( 136 ) 



Consider the straight line n + uik = as is shown on the lattice of (n, fc)-strings on fig. [|. This line is directed 
along the vector 

e(w) = uiei — &2 (137) 

orthogonal to the vector <j>, se e (p6| ). If charges of two BPS strings are both on the same side out of the line 
n + uk = then, according to ( |l3l| ) and (|8^), these strings do not interact and form marginally stable state. If 
instead they are from the opposite sides of this line, they attract each other and form a stable bound state. 

At u> = the straight line e(w) is directed along e 2 and as we increase u> moves anti-clockwise reaching the 
vector eo at u> — 1, see fig. |]. Moving it further, one reaches the vector ei at u> = oo. As the line directed along 
the vector e(w) hits any knot on the root lattice the BPS formula ( |87| ) 

T^ s = 2n£\n + ujk\ (138) 

gives zero for the BPS bound of corresponding string. Thus if within the two 2-7r/3 angles between the vectors 
— e 2 and ei (e 2 and — ei) on fig. ^ (see also "shadowed" regions at fig. |[ where this sector is shown explicitly) 
all stings were BPS saturated one would get infinitely many strings becoming tensionless. To avoid this we have 
to accept that all these strings are actually non-BPS states at real positive ui (cf. jn]). This is in accordance 
with our conclusion that BPS solutions exist only if n and k are both positive or negative, see (|S6"|). 

Thus, one gets the following picture for the (n, fc)-strings: ei- and e 2 -strings (— ei and — e 2 ), with charges 
proportional to the two closest roots on the SU(3) root lattice (see fig. |l|) are the lightest BPS stable " elementary" 
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CO- oo 




(0=0 



Figure 5: Different values of the parameter u>. Values < u> < 1 are restricted by straight lines uj — along the 
vector e2 and uj — I along eo, values lu > 1 are between the lines u> = 1 along eo and u> = oo, which is along 
the vector ei. 



strings. Both strings are from the same side of the line n + ajfc = 0, see fig. f|. Therefore all strings within two 
f angles between vectors ei and e2 (— ei and — e2) labeled by black circles at fig. [| exist as a BPS solutions of 
( |76|) but they are marginally unstable. Instead all strings labeled at fig. ||, by white circles can be considered as 
bound states of qe\ and — pe.2 strings (where q, p are integers and qp > 0), which are from different sides of the 
line n + ujk = 0. In this case the components are in attractive channel, hence, the composite strings are stable, 
but as we showed before they are not BPS saturated. 

Note as an additional check that (1 , 1) string is BPS marginally stable state. This is in complete agreement 
with our discussion at the end of sect. |4.2| where we have shown that this string should exist as a BPS solution 
of ([76]) in some region around u> = 1. Note also that e! and e 2 strings in fact do not become tensionless at 
uj — and uj — oo respectively as suggested by BPS formula (138). As we already mentioned "boundary" values 



uj = and uj = oo corresponds to strong coupling where one of the photons becomes massless and classical BPS 
mass formulas should be modified. 



Note also that there is a hidden symmetry of a W-boson spectrum ( 134 ) ^| 



M = fixed 

1 (139) 

z 

and let us check now that the spectrum of light particles in r — 2 vacuum also respects this symmetry. The 



eigenvalues of mass matrix are given by eqs. fl48|), (49), and in terms of variables (|135[) they get the form 




(mf)± = - 5 V»f \- + z±d^-l + z 2 } (140) 



10 This symmetry is very similar to a symmetry in Toda chain integrable system, which has already appeared in the context of 
SUSY gauge theories as an elegant way to formulate the Seiberg-Witten exact solution |36|. 
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Therefore we conclude that the low energy spectrum is also invariant under the symmetry (139). Now consider 
the spectrum of BPS (n, k) strings given by (138), which in terms of M and z (135) reads 



n 

- + kz 

z 



(141) 



We see that it is indeed invariant under transformation (13S) together with exchange in the string quantum 
numbers 

(n,k) <-> (k,n) (142) 

To conclude this section let us sum up the picture of the monopole confinement at r = 2 vacuum. The 
ei-monopole-antimonopole pair is confined by the ei-string, while e2-monopolc-antimonopole pair is confined 
by the e2-string. The eo-monopole-antimonopole pair is confined by the eo-string which is stable bound state 
of ei- and — e2-strings. We see that one gets three monopole-antimonopole meson states instead of one. As we 
have already explained, this "multiplicity" reflects the Abelian nature of confinement in Seiberg-Witten theory. 
For generic values of w the masses of all three mesons are different (they are determined by different values of 
string tension). 

Finally, let us note, that we predict presence of the non-BPS stable strings which can be considered as bound 
states of "elementary" BPS strings in attractive channel. These strings can form an "exotic" multi-monopole- 
multi-antimonopole mesons. The example of such string is (-1,2) string which is a stable non-BPS bound state 
of (— eo)- and e2-string. This string bound e2-monopole and eo-antimonopole with e2-antimonopole and eo- 
monopole to form an "exotic" meson. Presence of these "exotic" mesons also reflects the Abelian nature of 
confinement in Seiberg-Witten theory |37], 0. However, we do not expect such states to appear in theory with 
non- Abelian confinement. 



6 SU{2) gauge symmetry restoration in the theory with four and 
five flavors 

Finally, let us turn to the most physically interesting examples - the ST/(3) gauge theory with Nf = 4 and 
especially Nf — 5 flavors. As we already explained, the motivation to study such theories is that at certain 
submanifolds in their moduli space non-Abelian ST/ (2) gauge symmetry is restored and one can study what 
happens to confinement and flux tubes in this regime. The reason for restoration of non-Abelian SU (2) gauge 
symmetry in quantum theory is that at Nf = 4 and Nf = 5 the corresponding SU(2) subsectors are not 
asymptotically free and they remain to be non-Abelian in weak coupling regime at ttia S> A (cf. 0). 

Let us take a closer look at r = 2 vacua of these theories. The theory with Nf = 4 has six r = 2 vacua 
while the theory with Nf — 5 has ten r = 2 vacua, see (^lj). Now consider the point in the parameter space 
where masses of all flavors ttia become equal. At this point six (ten) vac ua o f the theory with Nf = 4 (5) 



flavors collide with the two-fold effect. First, as we already discussed in sect. |2.4| , the Higgs branch of dimension 
8 (Nf — 2) develops. It touches the Coulomb branch at the point 



a-3 



0, a 8 = -VEm, (143) 



where m = tiia is now common value of the quark masses, see (ff4|). 

Second, since (13 = the SU (2) gauge subgroup acting in the upper left corner of the adjoint field matrices 
is classically restored, see ( |3l| ) FJ. this restoration is also preserved in quantum theory. Namely, the SU(2)- 
subsector for Nf — 4 is conformal and the coupling constant does not run, being of the order of 

9 2 ~ 1 r m « 1 ( 144 ) 

log (m/ A) 

(with A = A-su(3)> the SU(3) theory with Nf — 4 is not conformal), frozen at the scale to where SC/(3) group 
is broken down to SU{2) x U(l). 

11 Strictly speaking this is true at energies much larger then quark VEV's which break gauge group completely at the scale of 
order of ^Jfim, see discussion below. 
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For Nf = 5 theory the £{7(2)-subsector has "zero charge" in the infrared limit. The coupling constant at 
large distances is of the order of 

9 1 

g : — n — ; , 9N < 1, (145) 

log(AmA/m 2 ) v ; 

frozen at lower scale of quark mass difference Am (note that Agy^) = m 2 / A.su(S) = rn 2 /A for Nf = 5). 

We consider here special submanifold of the Higgs branch which admits the BPS flux tubes (cf. [0, |Ts||). 
We call this submanifold origin of the Higgs branch. One point on this submanifold corresponds to non-zero 
u-component of the first flavor and non-zero d-component of the second flavor given by 

(inu 1 ) = (d 2 d 2 ) = 3/wn = | (146) 

while all other components are zero. This is a continuation of r = 2 vacuum considered in previous sections to 
the case of e qua l quark masses, see JlS|). Other points in the origin of Higgs branch are given by SU(Nf) flavor 
rotation of ( |l46[ ) . The dimension of the origin of the Higgs branch is 4(Afy- — 2). To see this note that VEV's of 
u and <i-quarks in r = 2 vacuum break SU(Nf) symmetry down to SU(Nf — 2). Thus the number of "broken" 
generators is dim SU(Nf) — dim SU(Nf — 2) = 4(7V/ — 1) and four phases are "eaten" by Higgs mechanism. 

Other points on the 8(Nf — 2) dimensional Higgs branch correspond to non-zero VEV's of massless moduli 
fields, and these points do not admit BPS strings. In particular, the ANO strings on Higgs branch were studied 
in ||39| , they correspond to limiting case of type I strings with the logarithmically thick tails associated with 
massless scalar fields. Moreover, the infinitely long strings do not exist |3^, and we do not discuss here 
strings in generic points on Higgs branch. 

Now let us focus on the strings arising in some particular vacuum at the origin of Higgs branch. By flavor 
rotation one can always put this vacuum to the form ( |146| ). To study the strings at this vacuum we can apply 
results of the previous section and take the limit lo — > 1 . 

Suppose we approach the value w = 1. Let us make a closer look at what happens to the lattice of strings 
from fig. |H It is easy to see that eo-string corresponds to winding around the {7(1) factor associated with the 



generator A 3 , see (149). This generator belongs to the restored £{7(2) subgroup. However, £{7(2) group does 
not admit flux tubes just because m(SU(2)) — 0. Now it is clear that e -string becomes unstable at u — > 1. In 
fact, the winding just shrinks to zero once the group manifold becomes 3-sphere instead of a circumference. 

What happens physically is that eo-string is broken by eo-monopole-antimonopole production which does 
not cost any energy in the limit of equal quark masses. To see this note that mass of the eo-monopole vanishes 
in this limit. Namely, the BPS eo-monopole mass in weak coupling regime has the form pj 

(147) 



and vanishes at C13 — > 0. 

Thus, we see that eo-string becomes unstable and disappears from the spectrum. Let us discuss now what 
happens to confined states when we tend to zero the differences of quark masses Atuab = tua — tub (assuming 
that they are all of the same order Attiab ~ Am). As we already explained at JJixh <C Am <C m one has 
wea k co upling and the mass of eo W-boson is of the order of Am, while the mass of eo-monopole is given 



by (147). When we put Am well below JJmi the mass of eo W-boson (together with the masses of photons) 
is determined by complete breaking of SU(2) gauge subgroup by the quark VEV's (instead of adjoint VEV 
(as) ~ Am), and it is frozen at the val ue of o rder of Jiim. As we reduce Am the eo-monopole becomes lighter 



than eo W-boson. This shows that eqs.(|144j), (145) are no longer valid and we enter into strong coupling regime. 
To understand this, note that presence of simultaneously light quarks (which become massless quark moduli at 
Am = 0) and light monopole means approaching the point of Argyres-Douglas type ||^, |[J, . Apparently 
as eo-monopole becomes lighter and lighter the eo-string becomes more and more unstable. 

As we already discussed, presence of N c — 1 different strings associated with each U(\) factor of SU(N C ) gauge 
group broken down to {/(l)^ " 1 leads to unwanted multiplicity in the hadron spectrum of Af = 2 QCD which 
we do not expect in ordinary QCD ||, and this reflects the {7(1) nature of the confinement in Seiberg-Witten 
theory. In particular, one has generically two "elementary" strings in £{7(3) gauge theory, namely ei-string 
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and e2-string in r = 2 vacuum. As we already explained in the last section this leads to the presence of three 
monopole-antimonopole meson states formed by ei, e.2 and eo strings. 

However, in the theory with Nf — 4 and Nf = 5 flavors with equal masses with the non-Abelian SU(2) 
subgroup of gauge symmetry restored at r = 2 vacuum eo-string becomes unstable. This has two-fold effect 
on the meson spectrum. First eo-monopole-antimonopole meson formed by eo-string with eo-monopole and 
antimonopole attached to its ends acquire large width. Eventually at Am — > it becomes unobservable as a 
resonance state and disappear from the spectrum. 

Second, ei and e2 monopole-antimonopole mesons formed by ei and e2 strings respectively become mixed 
with the transition amplitude of order of one because e! — e 2 = e . Hence, these two mesons form two splited 
states, more massive state acquires large width and also disappears from the spectrum. As a result we end up 
with only one set of monopole-antimonopole meson Regge trajectories. 

Hence, we see that there is no longer unwanted multiplicity in the hadron spectrum. Although confinement 
in this theory is due to the presence of Abelian ei-string (~ e2-string), the hadron spectrum has multiplicity 
which we expect in a theory with non-Abelian confinement. 

Note also that extra multi-monopole-muti-antimonopole "exotic" states formed by stable non-BPS strings 
(see sect. 5.3) also disappear in the equal mass limit in theories with Nf — 4 or Nf = 5. 

To conclude this section let us make a comment on so called semilocal strings (see [Q for a review). These 
arc string-like solutions which interpolate between BPS strings and two-dimensional sigma model instantons 
lifted to four dimensions. Instead of ordinary BPS strings these solutions have power fall-off at large r. The 
semilocal strings arise in the models with additional global symmetry. In the theories with Nf = 4 and Nf — 5 
considered in this section one has broken global SU(Nf) symmetry at the origin of Higgs branch which leads 
to presence of massless Goldstone scalars responsible for possible power behavior of string solution at large r. 
The presence of semilocal string manifests itself as a zero mode of the ordinary BPS string. 

Therefore it would be interesting to investigate how ordinary BPS strings with exponential fall-off at infinity 
studied in this paper turn into semilocal strings as we aproach the point of SU (2) gauge symmetry restoration. 
One comment here is that as it is mentioned in sect. |4 . 1| our strings are BPS saturated only in the leading 
order in fi/m when perturbation of superpotential (|]J) reduces to the FI term and Af = 2 supersymmetry is 
not broken. Already in the next to leading order in fx/m these strings turn into type I strings JlJ]. For type I 
strings the zero mode associated with the transverse size of the string (which is similar to the instanton size) 
becomes a positive mode. It means that the transverse size of the semilocal string although grows as compared 
to the size of the ordinary string (~ 1/^/im) is still bounded by l//! 3 / 4 ™ 1 / 4 , and the semilocal string does not 
grow infinitely thick. 



7 Conclusion 

In this paper we have studied the flux tubes in softly broken JV = 2 QCD arising in quark vacua of SU(N C ) 
gauge theory with Nf flavors at weak coupling, in particular focusing on the theory with SU (3) gauge group. 
These magnetic flux tubes are responsible for the confinement of monopoles. Of course, it would be more 
desirable to study confinement of quarks arising in dyon vacua at strong coupling via electric flux tubes. Still 
we believe that the two problems are directly related due to monodromies ||, |) and it is useful to start with 
simpler problem of monopole confinement in quark vacua, which one can keep at weak coupling. 

Generically £{7(3) gauge group is broken down to U(l) 2 , thus in each Af = 1 vacuum two "elementary" flux 
tubes arise. In r = 1 vacua (when only u-quark condense) these are magnetic u-string and electric e 2 -string. 
This hybrid phase arises since classically the gauge group is broken down only to SU(2) x U(l). However 
on quantum level the SU(2) factor is further broken down to another U(l) due to the Seiberg-Witten strong 
coupling mechanism. This results in the condensation of e2-monopole (dyon) and formation of electric e2-string. 
We have shown that strings in r = 1 vacua are the standard BPS ANO vortices in the limit of small adjoint 
masses \i. 

In r = 2 vacua with non-zero VEV's of u- and <i-quarks two "elementary" BPS strings are magnetic eo 
and e2-string (at the values < u> < 1 of VEV's ratio u> — (tferf 2 ) / (uiu 1 )). These strings are generalization of 
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standard ANO vortices, and they involve two gauge potentials interacting with two scalar fields satisfying BPS 
first order eqs. (|7^). We have demonstrated that bound states of multiple ei-string and multiple e2-string are 
marginally stable strings while bound states of multiple ei-string and multiple e2-antistring are stable non-BPS 
strings. 

Finally we have considered the theory with Nf — 4 and Nf = 5 in the limit of equal quark masses. In this 
limit SU (2) subgroup of the gauge group is not broken even on quantum level in r = 2 vacuum. We have shown 
that in the limit of equal quark masses eo-string becomes unstable and disappears from the spectrum. The 
monopole confinement is due to remaining ei-string which is not distinguishable in this limit from e2-string. This 
mechanism eliminates the unwanted multiplicity of the hadron spectrum generically present in Seiberg-Witten 
theory. The meson spectrum looks as expected in a theory with non-Abelian confinement. Namely, one gets the 
only set of ei-monopole-antimonopole meson Regge trajectories (or, equivalently, e2-monopole-antimonopole 
meson). 

Something peculiar happen to "baryons" in this limit. Of course we cannot really talk about baryons formed 
by monopoles since monopoles have zero baryon number. Still we can consider neutral baryon-like states formed 
by eo, e2 and §i = — ei monopoles connected by — ei and e2 strings in a chain-like configuration. Generically 
when SU(3) gauge group is broken down to U(l) 2 there are six different "baryon" states of this kind. However, in 
the limit of equal quark masses when SU (2) gauge subgroup is restored and e -string disappears these " baryons" 
reduce to the only meson formed by ei-monopole and ei-antimonopole. It is possible because monopoles have 
zero baryon number. 

This is not exactly what one expects in a theory with non-Abelian confinement. In a theory with non-Abelian 
confinement one expects to have one set of meson and one set of baryon Regge trajectories. In the theory at 
hand in the limit of equal quark masses " baryon" states disappear and we are left with only one set of meson 
trajectories. It would be interesting to study the mechanism of non-Abelian confinement via Abelian flux tubes 
suggested in this paper for the case of monopole condensation and quark confinement. In particular, this may 
allow us to find what happens with baryon multiplicity upon non-Abelian gauge subgroup restoration. 
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Appendix. SU(3) conventions and beta-functions 

The variables 4>% are associated with the basis vectors of the fundamental representation 3 (see fig. |l|) : fa = 4>fi { , 
i = 1,2, 3, where <p is a VEV vector in Cartan plane. Their relation with Cartesian co-ordinates <j> = (03, as) is 
given by 
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(148) 



i.e. ($) is decomposed over well-known diagonal Gell-Mann matrices 
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(149) 
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and 



A 8 = 
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^=(ai2 + a 2 ) = ^=(Mi + M 2 - 2m 3 ) 



(150) 



where we have used the following identification between the vectors in Cartan plane (see fig. |l|) and diagonal 
matrices 

/ 1 \ 
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V o o -i I 



= Mi - ^3 



(151) 



(152) 



(vectors (151) and (152) are normalized to a 2 = 2 and the factor | in (31), ( |32j ) is related to conventional 
normalization of the Gell-Mann matrices Tr(A a A&) = \8 a b). 

The W-boson masses ([71), (B8T) are proportional to 



acj) 



(153) 



are given then by the scalar products with the roots so that the (perturbative) prcpotential for pure gauge 
theory may be rewritten as 

OLCf) 

X 



E w>) 2 io g - 



aeA 



(154) 



+ 



where sum is taken over all positive roots A+. The quadratic part of (154) can be easily rewritten in terms of 
Cartesian co-ordinates or weights of the fundamental representation, using 



aeA + 



or 



^ 2 = ^E(^) 2 = ^rE 



1 JVe-l 
"E^ 



(155) 



(156) 



The matrix of coupling constants for the SU(3) gauge theory may be easily computed in terms of the root 
variables A\ = <f>\ — <fe = ati 2 4> an d A2 = 4>2 — <f>3 = OL 2 <p, where we used one of the simple roots a 2 and the 
"highest" root a 12 = a.i + a 2 (see fig. |l|) 
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and it can be rewritten in following co-ordinates 
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If all Ai ~ M -► oo 



61ogM 



1 
1 



0(1) 



ddiddj 

with the coefficient /3ym = 2N C = 6. 

One may also consider the contribution of the fundamental multiplets to the set of effective constants (|l5 
They come from the contribution of fundamental massive multiplets to the prepotential 
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whose second derivatives w.r.t. give rise to the contribution to (157) of the form 
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In following co-ordinates one gets 
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In the limit Ai ~ M — ► oo ( |163| ) gives rise to 
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with the coefficient giving necessary contribution into (3qcd = 2AT C — Afy. 
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